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ABSTRACT 

As  the  traditional  filtering  methods  may  cause  the 
problem  of  blur  when  applied  on  stripe-like  images, 
this  paper  raised  a  new  method  intending  to  solve  the 
problem,  and  provided  the  process  of  the  method 
development.  The  method,  based  on  genetic 
algorithms  as  well  as  two-dimensional  Fourier 
transform,  offered  a  new  perspective  to  denoise  an 
image  with  abundant  components  of  high  frequency, 
and  is  of  certain  significance. 

Keywords  :  2-D  Fourier  transform,  genetic 
algorithms,  autocorrelations 

I.  Operating  Environment 

This  method  is  developed  based  on  Matlab  2014b. 
Corresponding  or  updated  version  is  recommended 
when  verifying  to  ensure  that  the  program  runs 
normally.  As  the  algorithm  used  a  multi-cycle,  all 
three  versions  of  the  program  require  much  time.  For 
the  first  edition  of  the  program,  when  the  number  of 
individuals  is  only  10,  running  a  generation  requires 
up  to  10  minutes  or  so,  the  second  version  of  the 
program  also  takes  about  6  minutes  per  generation, 


modulus  of  spectrum  of  source  Image,  in  In 


while  the  third  edition  of  each  generation  about  2 
minutes. 

II.  Problem  background 

This  method  is  mainly  applied  to  the  repair  of  stripe¬ 
like  images  polluted  by  Gaussian  noise.  The  general 
idea  of  filtering  an  image  is  to  preserve  its  low 
frequency  component  meanwhile  remove  the  high 
frequency  component.  For  striped  images  (such  as 
zebra),  the  frequency  domain  of  high  frequency 
components  is  fairly  abundant.  The  problem  of  image 
blur  will  emerge  on  conditions  of  direct  filtering.  So  it 
is  necessary  to  adopt  a  new  method  for  filtering. 

In  the  preparation  of  each  version  of  the  program, 
ideas  have  undergone  major  changes,  such  as  the 
shape  of  the  filter.  So  parts  of  the  Notes  or  other 
statements  that  does  not  affect  the  results  may  be  left 
in  the  program. 

III.  Methods 

Firstly,  the  characteristics  of  the  original  image 
spectrum  and  the  noise  image  are  analyzed,  as  shown 
in  the  following  figure. 


modulus  of  spectrum  of  image  with  noise  of  20  dB,  in  In 


Figure  3.1  modulus  of  spectrum  with  or  without  noise 
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An  image  of  zebra  is  used  here.  From  the  nature  of  the 
Fourier  transform  we  know  that  the  spectrum  of  the 
noise  image  is  equal  to  the  simple  superposition  of  the 
noise  spectrum  and  the  image  spectrum,  i.e. 

F(A  +  n)  =  F(A)  +  F(n); 

Thus,  the  noise  spectrum  "submerges"  the  spectrum  of 
the  image.  Since  the  noise  spectrum  is  random  and 
unpredictable,  we  cannot  completely  restore  the 
image  spectrum  in  the  frequency  domain  when  the 
noise  spectrum  is  unknown,  but  only  have  the  result 
spectrum  approached  as  close  as  possible  to  the 
spectrum  of  the  original  image.  This  is  the  general 
idea  of  this  method.  From  the  above  graphs  we  can 
see  that  the  high-frequency  part  of  the  original  image 
is  obviously  subsided,  and  after  the  noise  is  added,  the 
high-frequency  region  is  obviously  uplifted  and 
relatively  flat.  In  the  low  frequency  part,  the  original 
image  spectrum  "terrain"  is  higher,  making  the 
"submerged"  effect  of  noise  not  significant.  Namely 
there  is  no  obvious  difference  in  the  trend  in  the  low 
frequency  part.  Therefore,  it  is  desirable  to  find  a 
filter  function  that  best  fits  the  image  through  a 
genetic  algorithm  so  that  it  can  properly  restore  the 
trend  of  the  spectrum. 

The  following  figure  shows  the  spectrum  of  the  image 
flattened: 


modulus  of  spectrum  of  source  image,  in  In 
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Figure  3.2  modulus  of  spectrum  of  source  image,  in  In 

In  order  to  make  the  algorithm  be  generic,  considering 
the  characteristics  of  the  image  spectrum,  the  filter  is 
set  to  a  series  of  concentric  ellipses.  In  order  to  obtain 
similar  attenuation  of  the  four  corners  of  the  high 


frequency  component,  it  is  also  required  a  gradually 
transition  from  ellipse  to  circle  as  the  major  axis 
increases.  The  parameters  of  the  filter,  i.e.,  the 
direction  of  the  ellipse,  the  initial  eccentricity,  the 
velocity  of  the  transition  to  the  circle,  the  trend  of  the 
attenuation  value,  are  used  as  the  genes  in  the  genetic 
algorithm,  and  the  optimal  value  is  given  by  the 
genetic  algorithm. 


example  filter 


Figure  3.3  example  filter 

In  particular,  taking  the  characteristics  of  abounding 
high  frequency  components  of  the  fringe  image  into 
account,  it  is  supposed  to  restore  the  high  frequency 
part  of  the  spectral  details  as  much  as  possible  under 
the  prerequisite  of  ensuring  that  the  image  spectrum 
trend  approaches  that  of  the  original  image.  For  the 
question  of  how  to  extract  the  details  of  high 
frequency  components  from  the  noise  image,  the 
method  of  obtaining  gradient  information  and 
denoising  with  threshold  is  adopted  in  this  method. 

The  gradient  method  is  applied  to  extract  the  edge 
information  of  the  image,  that  is,  the  high  frequency 
component  of  the  image.  In  this  method,  the  gradient 
of  the  image,  the  derivative  (the  difference)  of  the  two 
directions  x  and  y  ,  is  calculated  first,  giving  the 
direction  of  the  most  grayscale  change  at  each  point. 
The  directional  derivative  (difference)  is  obtained 
then  for  this  direction,  so  that  the  edge  information  of 
the  image  can  be  obtained.  When  the  image  is  added 
20dB  Noise,  the  image  edge  information  is  obtained  as 
follows: 
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edge  info  example 


Figure  3.4  example  of  edge  information 


In  order  to  restore  as  much  as  possible  the  noise-free 
high-frequency  components  from  the  information, 
considering  that  the  magnitude  of  "edge  information" 
formed  by  the  noise  tends  to  be  low,  a  threshold  is  set 
thusly,  ignoring  the  “edge  information”  where  the 
gray  value  is  less  than  the  threshold.  The  process 
effect  is  as  follows: 


edge  info  with  threshold  =  0.36 


Figure  3.5  filtered  edge  information  with  threshold  =  0.36 


It  can  be  seen  that  although  the  edge  of  the  image  is 
also  lost  in  some  degree,  it  still  retains  part  of  the  edge 
information.  With  this  part  of  the  edge  information 
(rather  than  nothing  at  all),  it  will  be  possible  to 
modify  the  high-frequency  spectrum,  and  thus  try  to 
maximize  the  restoration  of  the  details  of  the  image. 

In  the  process  of  designing  this  method,  there  are 
three  editions  of  program  formed.  The  core  idea  of  the 
three  editions  is  to  use  genetic  algorithm  to  achieve 
the  optimal  individual  of  filtering  effect  through 
several  iterations.  For  genetic  algorithms,  there  must 
be  indicators  that  can  be  used  to  quantify  the 
restoration  of  individuals.  And  the  main  difference 
among  the  three  editions  of  the  program  is  that  the 
indicator  is  selected  different.  The  versions  are 
described  here  separately. 

3.1  The  first  edition  of  the  program 

The  evaluation  indicator  of  the  first  edition  of  the 
program  is  the  sharpness  of  the  central  peak  of  the 
autocorrelation  function  of  the  whole  image.  The 
following  demonstrates  its  principles. 

First  the  concept  of  two-dimensional  autocorrelation 
function  is  explained.  For  m  x  n  order  matrix: 


Define  its  cyclic  shift  as: 
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Referred  to  as  Apq  .  Obviously,  the  following 
expression  is  true  for  p  and  q : 

\\p\<m. 

1M<»  ’ 


Then  the  elements  within  its  autocorrelation  matrix 
can  be  given  by  the  following  equation: 


R«  = 


A  A 


mn 


pq  ’ 


Where  the  symbol  "  •  "  represents  the  dot  product 
operation  for  the  matrix,  that  is,  the  sum  after  the 
multiplication  of  the  corresponding  term.  We  can  see 
that  the  above  equation  is  actually  an  average  value. 
For  convenience,  the  autocorrelation  matrix  is 
expressed  as: 

Ra=^(A  A  pq);  (Equation  3.1.1) 

For  white  noise,  that  is,  Gaussian  noise,  we  know  that 
its  autocorrelation  characteristics  are  very  sharp,  the 
value  of  the  center  peak  is  very  large,  whereas  the 
non-center  value  is  very  small,  as  shown  below: 

R  of  50  x  50  matrix  of  white  gauss  noise,  25dB 
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Figure  3.1.1  R  of  matrix  of  25dB  white  gauss  noise 

For  the  general  image,  the  autocorrelation  matrix  is 
relatively  gentle,  as  shown  in  the  following  figure: 


R  of  picture  of  zebra 
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Figure  3.1.2  R  of  the  example  image 

So  the  filtering  effect  can  be  evaluated  by  judging  the 
sharpness  of  the  matrix  R . 

The  next  step  is  to  prove  that  the  autocorrelation 
matrix  of  superposition  of  the  image  matrix  A  and 
noise  n  is  equivalent  in  shape  to  the  simple 
superposition  of  the  two  autocorrelation  matrices,  i.e. 

RA+„  =RA+Rn+C; 

Available  from  Equation  3.1.1: 

RA+n  =  E((A  +  n)  •  ( A  +  n)M ) 

=  E((A  +  n)  •  (Apq  +  npq )) 

=  E( A-A  +  nn  +  An  +n-A  ) 

v  pq  pq  pq  pq  ’ 

Note  that  the  image  matrix  and  noise  should  be 
completely  independent,  so  there  are: 

E(n  ■  A  pq)  =  E(npq  •  A)  =  E(n)  ■  E(A)  =  C; 

In  this  way,  the  autocorrelation  matrix  of  the  image 
matrix  with  noise  is  equal  to  the  autocorrelation 
matrix  of  the  original  image,  where  the 
autocorrelation  matrix  of  the  noise  is  superimposed  as 
well  as  the  whole  lifted  by  the  height  of  constant  C . 
Thus  the  shape  is  equal  to  the  simple  superposition  of 
both. 


In  order  to  describe  the  sharpness  of  the  spikes,  also 
take  the  gradient  method  to  obtain  the  maximum 
value  of  the  directional  derivative  in  the  direction  of 
gradient  in  the  matrix  R  .  And  then  the  ratio  of  it  to  the 
average  difference  degree  (i.e.,  standard  deviation)  of 
matrix  R  is  used  as  a  basis  for  measuring  the 
sharpness  of  the  spikes,  which  is: 


@  IJTSRD  |  Available  Online  @  www.ijtsrd.com  |  Volume -1  |  Issue  -  5  |  July- Aug  2017 


Page:  1145 


International  Journal  of  Trend  in  Scientific  Research  and  Development  (IJTSRD)  ISSN:  2456-6470 

max(R)  eliminating  the  need  for  time-consuming 

^  ~  a  ’  autocorrelation  operation. 


The  lower  the  ratio,  the  better  the  filtering  effect.  In 
the  genetic  algorithm  thusly  this  ratio  can  act  as  an 
evaluation  function,  ascending  order. 

3.2  The  second  edition  of  the  program 

Although  the  first  version  of  the  program  is 
theoretically  feasible,  but  in  practice  due  to  the 
autocorrelation  function  call  for  too  much  time, 
another  alternative  evaluation  indicators  is  required  as 
substitute.  The  second  edition  of  the  program 
therefore  takes  the  mean  of  power  spectrum  as  the 
evaluation  index. 

From  the  Wiener-Khinchin  theorem  we  know  that  the 
power  spectrum  and  the  autocorrelation  function  are  a 
pair  of  Fourier  transforms,  i.e. 

S  =  F(R); 

So  we  can  get  the  power  spectrum  of  the  noise  image: 

sAt.  =F(RA+„) 

=  F(Ra+R„+C) 

=  SA+S,+<? 

If  we  ignore  the  central  impact,  we  can  see  that  the 
power  spectrum  of  the  noise  image  is  the 
superposition  of  the  power  spectrum  of  the  original 
image  and  the  noise  image.  So  there  is: 

F(SA+n)  =E(F(  R 

A+n  )) 

=  E(SA)  +  E(Sn) 

And  because  of: 

SA+n  =|F(A  +  n)|2; 

There  must  be  a  non-negative  expected  value  of  the 
power  spectrum.  So  f(SA+n)  can  be  regarded  as  the 
evaluation  indicator.  The  lower  value  of  the  indicator 
indicates  more  ingredients  of  E( Sn)  filtered  in  the 
image,  i.e.  better  filtering  effect.  In  this  way,  only  the 
procedure  of  having  the  already  calculated  elements 
of  the  spectrum  squared  and  averaged  is  needed, 


In  fact,  this  version  of  the  program  has  an  erroneous 
principle.  Since  the  filter  is  generated  on  the  basis  of 
the  unknown  contamination  of  the  image,  the  filter 
generation  statement  cannot  distinguish  E( Sn)  versus 
E{ SA)  ,  but  only  blindly  pursue  the  minimum  of 
^(^A+n)  ,  leading  ingredients  of  E{ SA)  will  also 
suffer  a  lot  of  losses.  There  is  no  way  to  solve  this 
problem. 

3.3  The  third  edition  of  the  program 

The  first  two  versions  of  the  program  is  mainly 
aiming  at  finding  a  filter  for  the  image  through  the 
analysis  of  the  image,  in  order  to  achieve  optimal 
filtering.  Although  theoretically  feasible,  the  results  of 
the  operation  of  the  filtering  effect  is  very  poor,  and 
takes  much  time.  In  this  context,  the  third  version  of 
the  program's  aim  has  undergone  great  changes.  The 
third  edition  of  the  program  tries  to  find  a  filter  that  is 
relatively  effective  for  all  images  through  sample 
training.  In  order  to  make  the  filter  have  generality,  as 
well  as  to  compress  required  time  for  the  program,  a 
completely  concentric  model  is  adopted  instead  of  the 
ellipse-to-round  transition  filter  model.  Since  it  is  a 
sample  training,  the  evaluation  indicator  can  be 
performed  under  conditions  where  free  noise  images 
are  already  available.  Therefore,  the  mean  square 
error  of  the  filtered  image  and  the  original  image  is 
used  as  the  evaluation  indicator: 

m  n  0 

ZX(4-4) 

MSE  =  - ; 

mn 

IV.  Results  of  the  operation 

4.1  The  first  edition  of  the  program 

The  following  figure  is  the  result  of  the  first  version 
of  the  program  with  the  number  of  individual 
population  set  of  10,  generations  of  40,  mutation  rate 
of  0.1,  noise  of  25dB  .  The  program  requires  71.2 
minutes. 


@  IJTSRD  |  AvailableOnline@www.ijtsrd.coml  Volume  -  1  |  Issue  -  5  |  July- Aug  2017 


Page:  1146 


International  Journal  of  Trend  in  Scientific  Research  and  Development  (IJTSRD)  ISSN:  2456-6470 

O  Result  Monitor  —  □  X 
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Figure  4.1  result  of  the  first  edition  of  the  program 


4.2  The  second  edition  of  the  program 


The  following  figure  is  the  result  of  the  second  version  of  the  program  with  the  number  of  individual 
population  set  of  10,  generations  of  40,  mutation  rate  of  0.1,  noise  of  25dB.  The  program  requires  48.4  minutes. 


Figure  4.2  result  of  the  second  edition  of  the  program 

4.3  The  third  edition  of  the  program 

The  following  figure  is  the  result  of  the  third  version  of  the  program  with  the  number  of  individual  population 
set  of  10,  generations  of  40,  mutation  rate  of  0.1,  noise  of  25dB  .  The  program  requires  6  minutes. 
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Figure  4.3  result  of  the  third  edition  of  the  program 


V  Results  and  discussion 

As  can  be  seen  from  the  results  of  the  operation,  all  of 
the  three  editions  of  the  program  cannot  get  satisfying 
results  through  genetic  algorithm.  One  possible  reason 
is  that  because  of  the  limited  time,  parameter  of 
individual  number  is  set  too  small,  resulting  in  the 
case  where  genetic  algorithm  quickly  converge  to  the 
local  solution.  On  the  other  hand,  the  selection  of  the 
indicator  in  genetic  algorithm  selection  may  not  be 
scientific.  In  operation,  it  is  found  that  genetic 
algorithms  tend  to  screen  out  individuals  whose  edge 
information  is  completely  erased.  This  shows 
scientific  suspicion  of  the  measure  to  use  the 
unstandardized  gradient  edge  information  as  the 
replacement  of  high-frequency  component. 

Another  point  that  can  be  further  improved  is  the 
denoising  of  the  edge  information.  The  current  simple 
threshold  method  will  cause  a  lot  of  losses  of  edge 
information.  My  initial  idea  is  to  use  the  Hog 
algorithm,  that  is,  to  find  the  gradient  of  each  point  in 
the  image.  We  can  see  that  the  edge  of  the  noise  is 
often  expressed  as  a  scatter  without  continuity,  and 
the  edge  direction  of  the  real  image  is  within  a  certain 
range  consistent  with  the  direction  of  the  adjacent 
edge.  Therefore,  we  can  use  the  method  of  "moving 
hood"  in  median  filter  measure  to  analyze  the  gradient 
direction  in  a  small  area,  so  as  to  realize  the  filtering 
of  noise. 


Although  the  problem  of  Gaussian  denoising  of  high 
frequency  image  remained  unsolved,  this  paper 
presents  a  new  idea,  which  has  certain  significance. 

VI  Appendix:  Program  code 
6.1  The  first  edition  of  the  program 

%% 

%initializing 

clc; 

disp('GA  image  filtering  version  1.0'); 
disp(' '); 

disp('Initializing. . . ; 

clear; 

close  all; 

%% 

%reading  the  target  image 

disp('Loading  raw  image...'); 

source  =  double(imread('zebra/Zebra.png')); 

[XD,YD]  =  size(source); 
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max_r  =  sqrt((XD/2)A2  +(YD/2)A2); 

%max_S  =  XD*YD/4; 

%% 

%loading  parameters 
disp(' '); 

disp('Loading  parameters...'); 

generation  =  input('For  how  many  generations  are  you 
expecting?:  ');  %overall  generations  of  circulation  in 
GA 

while(generation  <=  1  ||  mod(generation,  1)  ~=  0) 

disp('The  parameter  is  supposed  to  be  a  positive 
integer  greater  than  T); 

generation  =  input('Please  re-input  the  parameter: 
end 

number  =  input('How  many  individuals  are  there  in 
the  flock?: ');  %number  of  individual 

while(number  <=  0  ||  mod(number,  2)  ~=  0) 

disp('The  parameter  is  supposed  to  be  a  positive 
even  number'); 

number  =  input('Please  re-input  the  parameter: '); 
end 

mutate  =  input('What  is  the  probability  of  mutation 
supposed  to  be?: ');  %probability  of  mutating 

while(mutate  <=  0  ||  mutate  >=  1) 

disp('The  parameter  is  supposed  to  be  located 
between  0  &  T); 

mutate  =  input('Please  re-input  the  parameter: '); 
end 

%% 

%adding  noise 
disp(' '); 

disp(' Adding  noise...'); 


dB  =  input('What  intensity  of  noise  are  you 
expecting?: '); 

noise  =  wgn(XD,  YD,  dB); 

test  =  double(source)  +noise; 

cal_t  =  fftshift(fft2(test)); 

angle_t  =  exp(li*angle(cal_t)); 

showt  =  log(l  +abs(cal_t)); 

%% 

%extracting  edge  information 

disp('Extracting  edge  information...'); 

[TestX,TestY]=gradient(test) ; 

bevel_t  =  sqrt(TestX.A2+TestY.A2); 

grad_t=TestX.  *(TestX./bevel_t)+TestY.  *(TestY./bev 
el_t); 

grad_t(isnan(grad_t))  =  0; 
clear  bevel  t  TestX  TestY; 

%% 

%GA  preparation 
disp(' '); 

disp('Initialization  accomplished,  preparing  for  the 
GA  circulation...'); 

filter  =  ones(size(source)); 

%the  contour  elipses  are  specified  by  the  following 
equation: 

%  (1  -c*cosA2(theta))*xA2  +(1  -c*sinA2(theta))*yA2 
+c*sin(2*theta)*x*y  =  constant 

%the  filter  function  is  as  following: 

%  h(x)  =  xA4  +para_l  *xA2 

%h(x)  ensures  value  at  the  center  of  the  filter  equals 
to  0,  which  will  contrbute  to  the  normalization 
procession. 

para  =  [0.7,  pi/3,  -1,  -10,  1,  0.36,  0.7]; 

%the  vector  of  parameter  is  sequenced  as  the 
following  order: 
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%  [para_e,  paratheta,  paravelocity,  para_2,  paral, 
pararange,  para_threshold] 

%where  para  range  indicates  the  range  of  the  filter's 
value,  parathreshold  restricts  the  threshold  for  the 
matrix  grad  to  be  set  to  zero,  as  well  as  para_velocity 
indicates  the  speed  from  a  elipse  to  transform  to  a 
circle 

limits  =  [0,  0.9999; 

0,  pi; 

-1,  0; 

-10, -2; 

0,  1; 

0,  1; 

0,  1]; 

limits_range  =  limits(:,  2)  -limits(:,  1); 

%the  limits  of  values  of  the  paras  are  as  follows: 

%  parae  paratheta  para_velocity  para_l 
para  range  para  threshold  para  coef 

%  [0,  1)  [0,  pi)  [-1,  0]  [-10,  -2]  [0,1]  [0,  1] 

[0,1] 

index  =  ones(l,  2*number); 

para  =  repmat(para',  1,  2*number); 

bestmark  =  nan(l,  generation); 

mark_range  =  [0,  1]; 

filter  =  repmat(filter,  1,  l,2*number); 

result  =  repmat(test,  1,1,  2*number); 

show  f  =  repmat(show_t,  1,  1,  2*number); 

R  =  xcorr2(test); 

R  =  R./std(R(:)); 

R_range  =  [min(min(R(:))),  max(max(R( :)))]; 

[DiffX,  DiffY]  =  gradient(R); 

bevelr  =  sqrt(DiffX.A2  +DiffY.A2); 

gradr  =  DiffX.  *(DiffX./bevel_r) 

+DiffY .  *  (DiffY  ,/bevel_r) ; 


grad_r(isnan(grad_r))  =  0; 

primitive_mark  =  max(grad_r(:))/std(grad_r(:)); 

R  =  repmat(R,  1,1,  2*number); 

mark  =  primitive_mark*ones(l,  2*number); 

clear  primitive_mark  Restimated 

disp('Preparation  accomplished,  Launching  the  GA 
circulation...'); 

disp(' '); 

figure('Menubar', 'none', 'Name', 'Result 
Monitor', 'NumberTitle', 'off); 

set(gcf,  'position',  [239,  344,  1083,  420]); 

subplot(l,  2,  1); 

handle_image  =  imshow(result(:, :,  1),  []); 
title('Filtering  Result'); 
subplot(2,  2,  2); 
handlewindow  =  gca; 

area(best_mark,  'EdgeColor',  [77/255,  144/255, 

21/255],  'FaceColor',  [49/255,  171/255,  118/255], 
'YDataSource',  'bestmark'); 

hold  on; 

plot(best_mark,  'x',  'Color',  [249/255,  71/255, 

71/255],  'YDataSource',  'best  mark',  'LineWidth',  1); 

hold  off; 

grid  on; 

axis([l,  generation,  mark_range(l),  mark_range(2)]); 
title('Trend  of  Best  Score  among  Generations'); 
xlabel('Generations'); 
ylabel('Marks'); 
subplot(2,  6,  10); 
colormap(gca,  parula); 

surf(fdter(:, :,  1),  'ZDataSource',  'fdter(:, :,  1)'); 
shading  interp; 
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alpha(0.7);  %% 


xlim([l,  XD]); 

%GA  circulation 

ylim([l,  YD]); 

disp('Starting  the  first  generation...'); 

zlim([0,  1]); 

tic; 

set(gca,  'xtick',  [],  'xticklabel',  []); 

for  t  =  1  :  generation 

set(gca,  'ytick',  [],  'yticklabel',  []); 

%generating  new  population 

camproj('perspective'); 

couple  =  randperm(number);  %free  mating 

title('Filter'); 

for  i  =  1  :  2  :  number 

subplot(2,  6,  11); 

for  node  =  1  :  size(para,  1) 

colormap(gca,  parula); 

%swap  genes 

mesh(show  f(:,  :,  1),  'ZDataSource',  'show  f(:,  :, 

i)’); 

if(round(rand)) 

para(node,  i  +number)  =  para(node,  i); 

xlim([l,  XD]); 

para(node,  i  +number  +1)  =  para(node,  i  +1); 

ylim([l,  YD]); 

else 

zlim( [min(show_t( : )) ,  max(show_t( :))] ) ; 

para(node,  i  +number)  =  para(node,  i  +1); 

set(gca,  'xtick',  [],  'xticklabel',  []); 

para(node,  i  +number  +1)  =  para(node,  i); 

set(gca,  'ytick',  [],  'yticklabel',  []); 

end 

camproj('perspective'); 

%mutation 

title('Spectrum'); 

if(rand  <  mutate) 

subplot(2,  6,  12); 

para(node,  i  -l-number)  =  limits(node,  1) 

colormap(gca,  parula); 

+limits_range(node)*rand; 

contour3(R(:, :,  1),  30, 'ZDataSource',  'R(:, :,  1)'); 

end 

box  off; 

if(rand  <  mutate) 

xlim([l,  size(R,  1)]); 

para(node,  i  +number  +1)  =  limits(node,  1) 

ylim([l,  size(R,  2)]); 

+limits_range(node)*rand; 

zlim(Rrange); 

end 

set(gca,  'xtick',  [],  'xticklabel',  []); 

end 

set(gca,  'ytick',  [],  'yticklabel',  []); 

clear  node 

camproj('perspective'); 

end 

title('Autocorrelation'); 

clear  i; 

drawnow; 

%characterizing  individuals 
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for  i  =  number  +1  :  2*number 
%filter  generating 
for  j  =  1  :  XD; 
for  k  =  1  :  YD; 

ellipse  =  (1  -para(l,  i)*cos(para(2,  i))A2)*(j  - 
XD/2)A2  +(1  -para(l,  i)*sin(para(2,  i))A2)*(k  - 
YD/2)A2  +para(l,  i)*sin(2*para(2,  i))*(]  -XD/2)*(k  - 
YD/2); 

ellipse  =  pi*ellipse/sqrt(l  -(para(l,  i)A2)); 
circle  =  pi*((j  -XD/2)A2  +(k  -YD/2)A2); 
r  =  sqrt(circle/pi); 
ratio  =  r/max_r; 

portion  =  para(3,  i)*ratioA2  -(para(3,  i)  +l)*ratio 

+i; 

element  =  (portion*  ellipse  +(1 

portion)*circle)/(max_r)A2/pi; 

filter(j,  k,  i)  =  elementA2  +para(4,  i)*element; 

end 

clear  k; 
end 
clear  j; 

filter(:,  i)  =  filter(:,  i)/max(max(abs(filter(:, 
i))))*para(5,  i)  +1; 

%profile  modifying 

gradtemp  =  gradt; 

grad_temp(grad_temp  <  (max(grad_temp(:)) 
+std(grad_temp(:)))*para(6,  i))  =  0; 

calh  =  fftshift(fft2(grad_temp)); 

show  h  =  para(7,  i)*log(l  +abs(cal_h)); 

%filtering 

show_f(:,  i)  =  show_t.*filter(:,  i)  +show_h.*(l  - 
filter(:,  i)); 

cal  f  =  (exp(show_f(:, i))  -l).*angle_t; 
result(:, i)  =  abs(ifft2(fftshift(cal_f))); 


%scoring 

R(:, i)  =  xcorr2(result(:, i)); 

Rtemp  =  R(:,  i); 

R(:, i)  =  R_temp./std(R_temp(:)); 

[DiffX,  DiffY]  =  gradient(R(:,  i)); 

bevelr  =  sqrt(DiffX.A2  +DiffY.A2); 

gradr  =  DiffX.  *(DiffX./bevel_r) 

+DiffY.  *  (DiffY./bevelr); 

grad_r(isnan(grad_r))  =  0; 

mark(i)  =  max(grad_r(:))/std(grad_r(:)); 
end 
clear  i; 

%natual  selecting 
[mark,  index]  =  sort(mark); 
para  =  para(:,  index); 
filter  =  filter(:, index); 

R  =  R(:, index); 
show  f  =  show_f(:, index); 
result  =  result(:, index); 
bestmark(t)  =  mark(l); 

%demonstrating  champion 
markrange  =  minmax(bestmark); 
mark_range(l)  =  mark_range(l)*0.99; 
mark_range(2)  =  mark_range(2)/0.99; 
set(handle_image,  'CData',  result(:, 1)); 
set(handle_window,  'YLim',  mark_range); 
reffeshdata; 
drawnow; 

%informing 

disp(['Generation  #',num2str(t),'  has  finished']); 
toe; 
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beep;  %% 


dispO  ’); 

%terminate  circle 
if(t>2) 

if(best_mark(t)  ==  best_mark(t  -1)  &&  best_mark(t 
-1)  ==  best_mark(t  -2)) 

break; 

end 

end 

end 

%% 

%display  result 

figure('Menubar', 'none', 'Name', 'Result 
comparasionVNumberTitle','off); 

set(gcf,  'position',  [239,  162,  1083,  601]); 

subplot(l,  2,  1); 

imshow(test,  [min(test(:)),  max(test(:))]); 
title('Input  Image'); 
subplot(l,  2,  2); 

imshow(result(:, :,  1),  [min(test(:)),  max(test(:))]); 
title('Filtering  Result'); 

%clear  ellipse  a  ratio  portion  element  r  circle  square; 

6.2  The  second  edition  of  the  program 

%% 

%initializing 

clc; 

disp('GA  image  filtering  version  2.0'); 
disp(' '); 

disp('Initializing. . . '); 

clear; 

close  all; 


%reading  the  target  image 

disp('Loading  raw  image...'); 

source  =  double(imread('zebra/Zebra.png')); 

[XD,YD]  =  size(source); 

maxr  =  sqrt((XD/2)A2  +(YD/2)A2); 

%% 

%adding  noise 
disp(' '); 

disp('Adding  noise...'); 

dB  =  input('What  intensity  of  noise  are  you 
expecting?: '); 

while(isempty(dB)) 

dB  =  input('Waiting  for  input...’); 

end 

noise  =  wgn(XD,  YD,  dB); 
test  =  double(source)  +noise; 
calt  =  fftshift(fft2(test)); 
anglet  =  exp(li*angle(cal_t)); 
show_t  =  log(l  +abs(cal_t)); 

%% 

%loading  parameters 
disp(' '); 

disp('Loading  parameters...'); 

generation  =  input('For  how  many  generations  are  you 
expecting?:  ');  %overall  generations  of  circulation  in 
GA 

while(generation  <=  1  ||  mod(generation,  1)  ~=  0) 

disp('The  parameter  is  supposed  to  be  a  positive 
integer  greater  than  1'); 

generation  =  input('Please  re-input  the  parameter: 
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end 

number  =  input('How  many  individuals  are  there  in 
the  flock?: ');  %number  of  individual 

while(number  <=  0  ||  mod(number,  2)  ~=  0) 

disp('The  parameter  is  supposed  to  be  a  positive 
even  number'); 

number  =  input('Please  re-input  the  parameter: '); 
end 

mutate  =  input('What  is  the  probability  of  mutation 
supposed  to  be?: ');  %probability  of  mutating 

while(mutate  <=  0  ||  mutate  >=  1) 

disp('The  parameter  is  supposed  to  be  located 
between  0  &  T); 

mutate  =  input('Please  re-input  the  parameter: '); 
end 

%% 

%extracting  edge  information 

disp('Extracting  edge  information...'); 

[TestX,TestY]=gradient(test) ; 

bevelt  =  sqrt(TestX.A2+TestY.A2); 

grad_t=TestX.*(TestX./bevel_t)+TestY.*(TestY./bev 

elj); 

grad_t(isnan(grad_t))  =  0; 
clear  bevelt  TestX  TestY; 

%% 

%GA  preparation 
disp(' '); 

disp('Initialization  accomplished,  preparing  for  the 
GA  circulation...'); 

filter  =  ones(size(source)); 

%the  contour  elipses  are  specified  by  the  following 
equation: 

%  (1  -c*cosA2(theta))*xA2  +(1  -c*sinA2(theta))*yA2 
+c*sin(2*theta)*x*y  =  constant 


%the  filter  function  is  as  following: 

%  h(x)  =  xA4  +para_l  *xA2 

%h(x)  ensures  value  at  the  center  of  the  filter  equals 
to  0,  which  will  contrbute  to  the  normalization 
procession. 

para  =  [0.7,  pi/3,  -1,  -10,  1,  0.36,  0.7]; 

%the  vector  of  parameter  is  sequenced  as  the 
following  order: 

%  [para_e,  para_theta,  para_velocity,  para_2,  para_l, 
pararange,  para_threshold] 

%where  para  range  indicates  the  range  of  the  filter's 
value,  para  threshold  restricts  the  threshold  for  the 
matrix  grad  to  be  set  to  zero,  as  well  as  para_velocity 
indicates  the  speed  from  a  elipse  to  transform  to  a 
circle 

limits  =  [0,  0.9999; 

0,  pi; 

-1,  0; 

-10, -2; 

0,  1; 

0,  1; 

0,  1]; 

limits_range  =  limits(:,  2)  -limits(:,  1); 

%the  limits  of  values  of  the  paras  are  as  follows: 

%  parae  para_theta  para_velocity  para_l 
para  range  para_threshold  para  coef 

%  [0,  1)  [0,  pi)  [-1,  0]  [-10,  -2]  [0,1]  [0,  1] 

[0,1] 

%para_coef  has  currently  been  removed 
index  =  ones(l,  2*number); 
para  =  repmat(para',  1,  2*number); 
best  mark  =  nan(l,  generation); 
mark_range  =  [0,  1]; 

filter  =  repmat(filter,  1,  l,2*number); 
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result  =  repmat(test,  1,1,  2*number); 

show  f  =  repmat(show_t,  1,1,  2*number); 

grad_show  =  repmat(grad_t,  1,1,  2*number); 

power  =  (abs(fftshift(fft2(test)))).A2; 

primitivemark  =  mean(power(:)); 

mark  =  primitive_mark*ones(l,  2*number); 

percentage  =  [0,  0]; 

clear  primitive  mark  Restimated 

disp('Preparation  accomplished,  Launching  the  GA 
circulation...'); 

disp(' '); 

figure('MenubarVnoneVNameVResult 
Monitor',  'NumberT  itle',  'off) ; 

set(gcf,  'position',  [237,  195,  1083,  650]); 

subplot(3,  5,  [1,  2,  6,  7]); 

handle_image  =  imshow(result(:, 1),  []); 

title('Filtering  Result'); 

subplot(3,  11,  23:32); 

handlewindow  =  gca; 

area(best_mark,  'EdgeColor',  [77/255,  144/255, 

21/255],  'FaceColor',  [49/255,  171/255,  118/255], 
'YDataSource',  'bestmark'); 

hold  on; 

plot(best_mark,  Y,  'Color',  [249/255,  71/255, 

71/255],  'YDataSource',  'best  mark',  'LineWidth',  1); 

hold  off; 

grid  on; 

axis([l,  generation,  mark_range(l),  mark_range(2)]); 
title('Trend  of  Best  Score  among  Generations'); 
xlabel('Generations'); 
ylabel('Marks'); 
subplot(3,  11,  33); 


area(percentage,  'EdgeColor',  'none',  'FaceColor', 
[49/255,  171/255,  118/255],  'YDataSource', 

'percentage'); 

set(gca,  'xtick',  [],  'xticklabel',  []); 
set(gca,  'ytick',  [],  'yticklabel',  []); 
ylim([0,  1]); 
title('Progress'); 
subplot(3,  5,  3); 
colormap(gca,  parula); 

surf(fdter(:, :,  1),  'ZDataSource',  'filter(:, :,  1)'); 

shading  interp; 

alpha(0.7); 

xlim([l,  XD]); 

ylim([l,  YD]); 

zlim([0,  1]); 

set(gca,  'xtick',  [],  'xticklabel',  []); 
set(gca,  'ytick',  [],  'yticklabel',  []); 
camproj('perspective'); 
title('Filter'); 
subplot(3,  5,  8); 
colormap(gca,  parula); 

mesh(show_f(:,  :,  1),  'ZDataSource',  'show_f(:,  :, 

i)’); 

xlim([l,  XD]); 
ylim([l,  YD]); 

zlim([min(show_t( :)),  max(sho  w_t( :))]); 
set(gca,  'xtick',  [],  'xticklabel',  []); 
set(gca,  'ytick',  [],  'yticklabel',  []); 
camproj('perspective'); 
title('Spectrum'); 
subplot(3,  5,  [4,  5,  9,  10]); 
handle  grad  =  imshow(grad_show(:, :,  1),  []); 
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title('Edge  Information');  end 


drawnow; 

clear  i; 

%% 

%characterizing  individuals 

%GA  circulation 

for  i  =  number  +1  :  2*number 

disp('Starting  the  first  generation...'); 

percentage  =  [1,  l]*(i  -number 

tic; 

1  )/number;refreshdata;drawnow; 

for  t  =  1  :  generation 

%filter  generating 

percentage  =  [0,  O];refreshdata;drawnow; 

for  j  =  1  :  XD; 

%generating  new  population 

for  k  =  1  :  YD; 

couple  =  randperm(number);  %free  mating 

ellipse  =  (1  -para(l,  i)*cos(para(2,  i))A2)*(j  - 
XD/2)A2  +(1  -para(l,  i)*sin(para(2,  i))A2)*(k  - 

for  i  =  1  :  2  :  number 

YD/2)A2  +para(l,  i)*sin(2*para(2,  i))*(j  -XD/2)*(k  - 
YD/2); 

for  node  =  1  :  size(para,  1) 

ellipse  =  pi*ellipse/sqrt(l  -(para(l,  i)A2)); 

%swap  genes 

circle  =  pi*((j  -XD/2)A2  +(k  -YD/2)A2); 

if(round(rand)) 

%square  =  ((max(abs(j  -XD/2),  abs(k  - 

para(node,  i  +number)  =  para(node,  i); 

YD/2)))A2); 

para(node,  i  +number  +1)  =  para(node,  i  +1); 

r  =  sqrt(circle/pi); 

else 

ratio  =  r/max_r; 

para(node,  i  +number)  =  para(node,  i  +1); 

%ratio  =  sqrt(square/max_S); 

para(node,  i  +number  +1)  =  para(node,  i); 

portion  =  para(3,  i)*ratioA2  -(para(3,  i)  +l)*ratio 

+i; 

end 

%mutation 

element  =  (portion*ellipse  +(1 

portion)*circle)/(max_r)A2/pi; 

if(rand  <  mutate) 

%element  =  (portion*  elipse/maxelipse  +(1  - 

portion)  *  square/ maxS) ; 

para(node,  1  -i-number)  =  limits(node,  1) 

+limits_range(node)*rand;  filter(j,  k,  i)  =  elementA2  +para(4,  i)*element; 

end  end 

if(rand  <  mutate)  clear  k; 

para(node,  i  +number  +1)  =  limits(node,  1)  end 
+limits_range(node)*rand; 


end 

clear  j; 

end 

filter(:,  :,  i)  =  filter(:,  :,  i)/max(max(abs(filter(:,  :, 
i))))*para(5,  i)  +1; 

clear  node 

%profile  modifying 
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grad_temp  =  gradt;  reffeshdata; 

grad  temp(grad_temp  <  (max(grad_temp(:))  drawnow; 


+std(grad_temp(:)))*para(6,  i))  =  0; 
grad_show(:, i)  =  grad_temp; 
calh  =  fftshift(fft2(grad_temp)); 
show_h  =  log(l  +abs(cal_h)); 

%filtering 

show_f(:,  i)  =  show_t.*filter(:,  i)  +show_h.*(  1  - 
filterO, i)); 

cal  f  =  (exp(show_f(:, i))  -l).*angle_t; 
result(:, i)  =  abs(ifft2(fftshift(cal_f))); 

%scoring 

power  =  (abs(cal_f)).A2; 
mark(i)  =  mean(power(:)); 
end 
clear  i; 

%natual  selecting 
[mark,  index]  =  sort(mark); 
para  =  para(:,  index); 
filter  =  filter(:, index); 
gradshow  =  grad_show(:, index); 
show_f  =  show_f(:, index); 
result  =  result(:, index); 
bestmark(t)  =  mark(l); 

%demonstrating  champion 
mark_range  =  minmax(best_mark); 
mark_range(l)  =  mark_range(l)*0.99; 
mark_range(2)  =  mark_range(2)/0.99; 
set(handle_image,  'CData',  result(:, 1)); 
set(handle_grad,  'CData',  grad_show(:, :,  1)); 
set(handle_window,  'YLim',  markrange); 


%informing 

disp(['Generation  #',num2str(t),'  has  finished']); 

toe; 

beep; 

disp(' '); 

%terminate  circle 
if(t>2) 

if(best_mark(t)  ==  best_mark(t  -1)  &&  best_mark(t 
-1)  ==  best_mark(t  -2)) 

break; 

end 

end 

end 

%% 

%display  result 

disp('Circulation  terminated'); 

figure('Menubar', 'none', 'Name', 'Result 
comparasion','NumberTitle','off); 

set(gcf,  'position',  [237,  195,  1083,  650]); 

subplot(l,  2,  1); 

imshow(test,  [min(test(:)),  max(test(:))]); 
title('Input  Image'); 
subplot(l,  2,  2); 

imshow(result(:, :,  1),  [min(test(:)),  max(test(:))]); 
title('Filtering  Result'); 

%clear  ellipse  a  ratio  portion  element  r  circle  square; 

6.3  The  third  edition  of  the  program 

%% 

%initializing 
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clc;  %% 

disp('GA  image  filtering  version  3.0');  %loading  parameters 


disp(' '); 


disp(' '); 


disp('Initializing. . . ; 

clear; 

close  all; 


disp('Loading  parameters...'); 

generation  =  input('For  how  many  generations  are  you 
expecting?:  ');  %overall  generations  of  circulation  in 
GA 


%% 

%reading  the  target  image 

disp('Loading  raw  image...'); 

source  =  double(imread(’zebra/Zebra.png')); 

[XD,YD]  =  size(source); 
disp('Loading  training  sample...'); 

=  double(imread('zebra/ZebraG 


test(:,  :,  1)  = 

20  22. 1947 .png')); 

test(:,  :,  2)  = 

30_18.8633.png')); 

test(:,  :,  3)  = 

40_16.5752.png')); 

test(:,  4)  = 

50_14.8842.png')); 

test(:,  :,  5)  = 

60_13.5614.png')); 


double(imread('zebra/ZebraG 

double(imread('zebra/ZebraG 

double(imread('zebra/ZebraG 

double(imread('zebra/ZebraG 


sample  =  size(test,  3) 
calt  =  test; 
anglet  =  test; 
showt  =  test; 
for  i  =  1  :  sample; 

cal_t(:, :,  i)  =  fftshift(fft2(test(:, :,  i))); 
angle_t(:, :,  i)  =  exp(li*angle(cal_t(:, :,  i))); 
show_t(:, :,  i)  =  log(l  +abs(cal_t(:, :,  i))); 
end 
clear  i; 


while(generation  <=  1  ||  mod(generation,  1)  ~=  0) 

disp('The  parameter  is  supposed  to  be  a  positive 
integer  greater  than  T); 

generation  =  input('Please  re-input  the  parameter: 
end 

number  =  input('How  many  individuals  are  there  in 
the  flock?: ');  %number  of  individual 

while(number  <=  0  ||  mod(number,  2)  ~=  0) 

disp('The  parameter  is  supposed  to  be  a  positive 
even  number'); 

number  =  input('Please  re-input  the  parameter: '); 
end 

mutate  =  input('What  is  the  probability  of  mutation 
supposed  to  be?: ');  %probability  of  mutating 

while(mutate  <=  0  ||  mutate  >=  1) 

disp('The  parameter  is  supposed  to  be  located 
between  0  &  T); 

mutate  =  input('Please  re-input  the  parameter: '); 
end 

%% 

%extracting  edge  information 
disp(' '); 

disp('Extracting  edge  information...’); 
grad_t  =  test; 
for  i  =  1  :  sample; 

[TestX,TestY]=gradient(test(:, :,  i)); 
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bevel  t  =  sqrt(TestX.A2+TestY.A2);  index  =  ones(l,  2*number); 


grad_t(:, 

i)=TestX.*(TestX./bevel_t)+TestY.  *(TestY./bevel_t); 
end 

grad_t(isnan(grad_t))  =  0; 
clear  bevelt  TestX  TestY  i; 

%% 


para  =  repmat(para',  1,  2*number); 
best  mark  =  nan(l,  generation); 
mark_temp  =  zeros(l,  sample); 
mark_range  =  [0,  1]; 
filter  =  repmat(filter,  1,  l,2*number); 


%GA  preparation 
disp('  ’); 

disp('Initialization  accomplished,  preparing  for  the 
GA  circulation...'); 

filter  =  ones(size(source)); 

%the  filter  function  is  as  following: 

%  h(x)  =  xA4  +para_l  *xA2 

%h(x)  ensures  value  at  the  center  of  the  filter  equals 
to  0,  which  will  contrbute  to  the  normalization 
procession. 

para  =  [-10,  0.4,  0.36,  0.42]; 

%the  vector  of  parameter  is  sequenced  as  the 
following  order: 

%  [para_velocity,  para_l,  pararange, 

parathreshold,  para_coef] 

%where  para  range  indicates  the  range  of  the  filter's 
value,  para  threshold  restricts  the  threshold  for  the 
matrix  grad  to  be  set  to  zero,  as  well  as  para_velocity 
indicates  the  speed  from  a  elipse  to  transform  to  a 
circle 

limits  =  [-10,  0; 

0,  1; 

0,  1; 

0,  1]; 

limits_range  =  limits(:,  2)  -limits(:,  1); 

%the  limits  of  values  of  the  paras  are  as  follows: 

%  para_l  para  range  para  threshold  paracoef 

%  [-10,-2]  [0,1]  [0,1]  [0,1] 


result  =  repmat(test,  1,  1,  1,  2*number); 
show_f  =  repmat(show_t,  1,  1,  1,  2*number); 
grad_show  =  repmat(grad_t,  1,  1,  1,  2*number); 
show  h  =  zeros(XD,  YD,  sample); 
for  j  =  1  :  sample 

diff_temp  =  result(:, :,  j,  1)  -  source; 
mark_temp(j)  =  std(diff_temp(:)); 
end 
clear  j; 

mark  =  mean(mark_temp)*ones(l,  2*number); 
percentage  =  [0,  0]; 

disp('Preparation  accomplished,  Launching  the  GA 
circulation...'); 

disp(' '); 

show  result  =  zeros(XD,  YD); 
show  grad  =  zeros(XD,  YD); 
f  show  =  zeros(XD,  YD); 
for  j  =  1  :  sample 

show  result  =  show_result  +  result(:, :,  j,  1); 
show_grad  =  show  grad  +  grad_show(:, :,  j,  1); 
f_show  =  f_show  +  show_f(:, :,  j,  1); 
end 
clear  j; 

show_result  =  show_result./sample; 
show_grad  =  showgrad./sample; 
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f_show  =  f_show./sample; 

figure('Menubar', 'none', 'Name', 'Result 
Monitor',  'NumberT  itle',  'off) ; 

set(gcf,  'position',  [237,  195,  1083,  650]); 

subplot(3,  5,  [1,  2,  6,  7]); 

handleimage  =  imshow(show_result,  []); 

title('Filtering  Result'); 

subplot(3,  11,23:32); 

handlewindow  =  gca; 

area(best_mark,  'EdgeColor',  [77/255,  144/255, 

21/255],  'FaceColor',  [49/255,  171/255,  118/255], 
'YDataSource',  'bestmark'); 

hold  on; 

plot(best_mark,  'x',  'Color',  [249/255,  71/255, 

71/255],  'YDataSource',  'best  mark',  'LineWidth',  1); 

hold  off; 

grid  on; 

axis([l,  generation,  mark_range(l),  mark_range(2)]); 
title('Trend  of  Best  Score  among  Generations'); 
xlabel('Generations'); 
ylabel('Marks'); 
subplot(3,  11,  33); 

area(percentage,  'EdgeColor',  'none',  'FaceColor', 
[49/255,  171/255,  118/255],  'YDataSource', 

'percentage'); 

set(gca,  'xtick',  [],  'xticklabel',  []); 
set(gca,  'ytick',  [],  'yticklabel',  []); 
ylim([0,  1]); 
title('Progress'); 
subplot(3,  5,  3); 
colormap(gca,  parula); 

%contour3(filter(:, :,  1),  30,  'ZDataSource',  'filter(:,  :, 

i)'); 


%box  off; 

surf(filter(:, :,  1),  'ZDataSource',  'filter(:, :,  1)'); 

shading  interp; 

alpha(0.7); 

xlim([l,  XD]); 

ylim([l,  YD]); 

zlim([0,  1]); 

set(gca,  'xtick',  [],  'xticklabel',  []); 
set(gca,  'ytick',  [],  'yticklabel',  []); 
camproj('perspective'); 
title('Filter'); 
subplot(3,  5,  8); 
colormap(gca,  parula); 
mesh(f_show,  'ZDataSource',  'f  show'); 
xlim([l,  XD]); 
ylim([l,  YD]); 

zlim([min(show_t( :)),  max(sho  w_t( :))]); 
set(gca,  'xtick',  [],  'xticklabel',  []); 
set(gca,  'ytick',  [],  'yticklabel',  []); 
camproj('perspective'); 
title('Spectrum'); 
subplot(3,  5,  [4,  5,  9,  10]); 
handle  grad  =  imshow(show_grad,  []); 
title('Edge  Information'); 
drawnow; 

%% 

%GA  circulation 

disp('Starting  the  first  generation...'); 
tic; 

for  t  =  1  :  generation 
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percentage  =  [0,  O];refreshdata;drawnow; 
%generating  new  population 
couple  =  randperm(number);  %free  mating 
for  i  =  1  :  2  :  number 
for  node  =  1  :  size(para,  1) 

%swap  genes 
if(round(rand)) 

para(node,  i  -i-number)  =  para(node,  i); 
para(node,  i  -i-number  +1)  =  para(node,  i  +1); 
else 

para(node,  i  -i-number)  =  para(node,  i  +1); 
para(node,  i  +number  +1)  =  para(node,  i); 
end 

%mutation 
if(rand  <  mutate) 

para(node,  i  +number)  =  limits(node,  1) 
+limits_range(node)*rand; 

end 

if(rand  <  mutate) 

para(node,  i  -l-number  +1)  =  limits(node,  1) 
+limits_range(node)*rand; 

end 

end 

end 

%characterizing  individuals 

for  i  =  number  +1  :  2*number 

percentage  =  [1,  l]*(i 

number)/number;refreshdata;drawnow; 

%filter  generating 

for  j  =  1  :  XD; 

for  k  =  1  :  YD; 

circle  =  pi*((j  -XD/2)A2  +(k  -YD/2)A2); 


filter(j,  k,  i)  =  -circleA2  +para(l,  i)*  circle; 
end 
end 

filter(:,  i)  =  filter(:,  i)/max(max(abs(filter(:, 
i))))*para(2,  i)  +1; 

%profile  modifying 

for  j  =  1  :  sample 

gradtemp  =  grad_t(:, j); 

grad_temp(grad_temp  <  (max(grad_temp(:)) 
+std(grad_temp(:)))*para(3,  i))  =  0; 

grad_show(:, j,  i)  =  grad  temp; 

calh  =  fftshift(fft2(grad_temp)); 

show_h(:, j)  =  log(l  +abs(cal_h)); 

end 

%filtering 
for  j  =  1  :  sample 

show_f(:,  j,  i)  =  show_t(:,  j).*filter(:,  i) 
+show_h(:, j).*(l  -filter(:, i)).*para(4,  i); 

cal_f  =  (exp(show_f(:, j,  i))  -l).*angle_t(:, j); 

result(:, j,  i)  =  abs(ifft2(fftshift(cal_f))); 

end 

%scoring 
for  j  =  1  :  sample 

diff_temp  =  result(:, j,  i)  -  source; 
mark_temp(j)  =  std(diff_temp(:)); 
end 

mark(i)  =  mean(mark_temp); 
end 

%natual  selecting 
[mark,  index]  =  sort(mark); 
para  =  para(:,  index); 
filter  =  filter(:, index); 
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grad_show  =  grad_show(:, index);  if(best_mark(t)  ==  best_mark(t  -1)  &&  best_mark(t 


show_f  =  show_f(:, index); 

-1)  ==  best_mark(t  -2)) 

result  =  result(:, :, :,  index); 

break; 

bestrnark(t)  =  mark(l); 

end 

%demonstrating  champion 

end 

mark_range  =  minmax(bestmark); 

end 

mark_range(l)  =  mark_range(l)*0.99; 

%% 

mark_range(2)  =  mark_range(2)/0.99; 

%display  result 

show  result  =  zeros(XD,  YD); 

disp('Circulation  terminated'); 

show  grad  =  zeros(XD,  YD); 

figure('Menubar','none', 'Name', 'Result 
comparasion','NumberTitle','off); 

for  j  =  1  :  sample 

set(gcf,  'position',  [237,  195,  1083,  650]); 

show  result  =  show_result  +  result(:, :,  j,  1); 

subplot(l,  2,  1); 

show_grad  =  show_grad  +  grad_show(:, :,  j,  1); 

imshow(source,  []); 

f_show  =  f_show  +  show_f(:, :,  j,  1); 

title('Target  Image'); 

end 

subplot(l,  2,  2); 

showresult  =  showresult./sample; 

imshow(show_result,  []); 

show_grad  =  showgrad./sample; 

title('Filtering  Result'); 

f_show  =  f_show./sample; 

%clear  ellipse  a  ratio  portion  element  r  circle  square; 

set(handle_image,  'CData',  showresult); 

set(handle_grad,  'CData',  show_grad); 

set(handle_window,  'YLim',  markrange); 

reffeshdata; 

drawnow; 

%informing 

disp(['Generation  #',num2str(t),'  has  finished']); 

toe; 

beep; 

disp(' '); 

%terminate  circle 

if(t>2) 
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